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We show that in the presence of ferromagnetic electronic reservoirs and spin-dependent tunnel 
couplings, molecular vibrations in nonmagnetic single molecular transistors induce an effective in¬ 
tramolecular exchange magnetic field. It generates a finite spin-accumulation and -precession for 
the electrons confined on the molecular bridge and occurs under (non)equilibrium conditions. The 
effective exchange magnetic field is calculated here to lowest order in the tunnel coupling for a 
nonequilibrium transport setup. Coulomb interaction between electrons is taken into account as 
well as a finite electron-phonon coupling. We show that for realistic physical parameters, an effec¬ 
tive spin-phonon coupling emerges. It is induced by quantum many-body interactions, which are 
either electron-phonon or Coulomb-like. We investigate the precession and accumulation of the con¬ 
fined spins as function of bias- and gate-voltages as well as their dependence on the angle enclosed 
by the magnetizations between the left and right reservoir. 
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I. INTRODUCTION 

Understanding the interplay between Coulomb inter¬ 
action and the coupling of the electronic charge to vi¬ 
brational degrees of freedom has motivated intense re¬ 
search in several areas of condensed matter physics in 
the recent years. Transistors in future nanoelectronic 
applications will ideally work with only a few physical 
carriers of information. Inherently, the dimension of the 
devices are scaled down, such that response times shorten 
and fast switching between the two states of a transistor 
(or a quantum dot) becomes possible. For few electron 
quantum dots, the relevant physical properties, such as 
the number of confined electrons, effective factors, and 
spin-orbit coupling strengths, are coherently controllable 
by tuning respective gate voltages in experiments, see 
Refs. 1 and 2 and references therein. Suitable physical 
setups include semiconducting heterostructures as well as 
carbon nanotubes^ or gated nanowires^ in order to con¬ 
fine a small number of electrons. Mechanical degrees of 
freedom are inherent to the particular geometry of these 
devices. 

In molecular transistors (or nano junctions), the elec¬ 
tronic and/or spin states of a molecule are used in a 
controlled way to realize logical operations'’^. Besides 
their electronic properties, molecules inherently possess 
significant internal vibrational degrees of freedom. Char¬ 
acteristic frequencies of internal mechanical vibrational 
modes of a molecule are used to witness a working de¬ 
vice in a transport setup^’^. In general, nano mechanical 
systems designed by means of a clamped nanobeam, ex¬ 
hibit vibrational degrees of freedom as well^. By properly 
adjusting electrostatic gates, nanomechanical quantum 
dots with a vibrational degree of freedom are created. 
Depending on the particular situation, the energy scales 
of the electronic energies and the vibrational motion may 
compete and may generate rich physical phenomena. For 
example, electronic energies can be significantly shifted 


due to a strong electron-vibrational coupling. This gives 
rise to the Franck-Condon blockade as observed in car¬ 
bon nanotubes^^, as well as a strong coupling between 
the tunneling of electrons and the mechanical motion^^. 

In addition to controlling charge states, the electronic 
spin degree of freedom has been addressed as well in elec¬ 
tronic devices^^. In quantum dots, electron spin reso¬ 
nance measurements allow to determine the particular 
state of the confined particle. Also in clean nanotubes, 
spin-orbit coupling effects serve to further split the single¬ 
particle energy spectrum in the absence of an external 
magnetic field. There, spin and valley degrees of freedom 
may be used to define appropriate qubits^’Since the 
electron spin is a natural two state system, control and 
coherent manipulation of a system with a few confined 
spins is desirable for applications in spintronics and for 
encoding quantum information. 

As a generic theoretical model for all the mentioned 
scenarios, the Anderson-Holstein model has been estab¬ 
lished in the recent years. A quantum dot or a molecular 
nano junction with a few electronic orbitals and a finite 
electron-phonon coupling is tunnel-coupled to large reser¬ 
voirs of free electrons. The interplay between the elec¬ 
tronic states and the mechanical motion in Anderson- 
Holstein setups has been in the focus recently. Franck- 
Condon blockade effects together with large Fano fac¬ 
tors have been studiedas well as the transport prop¬ 
erties of the system in the presence of a finite tempera¬ 
ture bias^^’^^. Both cases were analyzed in the weak-to- 
intermediate coupling regime. Thereby, tunneling pro¬ 
cesses up to second order in the dot-lead hybridization 
have been taken into account. For slow phonons (i.e., 
small phonon frequencies), the adiabatic phonon regime 
has been addressed^^’^^. Also in this regime, Jovchev et 
al}^ have performed a scattering state renormalization 
group study in order to compute the I{V) characteris¬ 
tics and to determine the spectral function under the 
condition of thermal equilibrium. An advanced pertur- 
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bative treatment of the electron-phonon coupling under 
nonequilibrium conditions has also been worked out^^. 
Likewise, Monte Carlo simulations have been performed 
in the transient regime^^. Restricting the Hilbert space 
to the subspace of a singly occupied electronic device, 
Flensberg et have used an incoherent rate equa¬ 

tion approach together with the assumption that the 
phonon is in thermal equilibrium to reproduce the main 
experimental features on a qualitative level^^. The in¬ 
fluence on transport of a nonadiabatic phonon has been 
studied in Ref. 15 and 25. In addition, the crossover 
between the regime of the adiabatic and the weak tun¬ 
neling has been addressed by the iterative summation of 
nonequilibrium quantum path integrals in Ref. 26, sug¬ 
gesting that the Franck-Condon blockade arises also at 
low temperatures. 

Further interesting effects come into play when quan¬ 
tum dots are contacted by ferromagnetic leads. Then, 
the reservoirs can induce effective exchange magnetic 
fields for the electrons in the quantum dot. This has 
become known as the spin-valve effect. Here, the elec¬ 
tronic spin quantum number is no longer a good quan¬ 
tum number, and the confined spins are allowed to pre- 
cess coherently^^^^^. Ferromagnets may be character¬ 
ized by their difference in the density of states of the 
majority and the minority spins. When two reservoirs 
are present (say, the left and the right one) with dif¬ 
ferent magnetizations, the tunneling magneto-resistance 
(TMR) depends sensitively on the angle enclosed by the 
directions of the two magnetizations. Certainly, this 
opens up the possibility to manipulate spin states and 
initialize or read out the particular states. A proper the¬ 
oretical description has to include a description of the 
quantum coherences appearing in the system. In the 
present context, quantum coherences arise due to the 
coherent spin evolution. A detailed study of the coher¬ 
ent time evolution of the electron spin in the presence 
of an electron-phonon coupling has not been reported 
so far in the literature. The influence of the electron- 
phonon coupling in quantum transport setups with fer¬ 
romagnetic leads has been studied by means of nonequi- 
librium Green’s function for the two limiting cases of van¬ 
ishing Coulomb interaction ([/ = 0) and infinitely large 
Coulomb interaction {U oo) in Refs. 30 and 31, re¬ 
spectively. Also the case of a single coupled ferromagnet 
used as a switch for the current polarization has been 
investigated recently^^. A Hartree-Fock-like approxima¬ 
tion is used to close the set of equations for the causal 
electronic Green’s functions, yielding the tunneling cur¬ 
rent and the tunneling magneto-resistance^^. We note 
that within these approximations, a quantum coherent 
spin evolution on the quantum dot is not obtained due 
to the absence of inhomogeneous magnetic fields and/or 
approximations regarding the Coulomb interaction be¬ 
tween electrons. In the context of nanocooling a quantum 
dot by spin-polarized nonequilibrium currents, a simi¬ 
lar setup as suggested here has been investigated^^’^^. 
There, the focus has been on the dynamical cooling and 



FIG. 1. (a) Sketch of the Anderson-Holstein spin-valve geom¬ 

etry. The spinful electronic level with energy Eo- is coupled to 
a single mechanical harmonic mode with frequency Q with a 
coupling constant A. The electrons are able to tunnel from the 
spin-polarized leads to the non-magnetic device due to finite 
tunneling couplings The left and right ferromagnets 

have different magnetizations, which are characterized by an 
angle = L/R. (b) The underlying coordinate system 

where 4> is the relative angle between the two directions Wl/r 
of the left and the right magnetization. 


heating of the phonon mode in presence of a magnetic 
molecular junction. 


In this work, we consider a nonmagnetic molecular 
bridge with weak coupling to ferromagnetic leads. We 
include, besides the electronic degrees of freedom, also a 
coupling of the electrons to a single mode of a molecular 
vibration. Furthermore, we restrict our analysis to the se¬ 
quential tunneling regime, which we explore by means of 
a real-time diagrammatic approach for the reduced den¬ 
sity matrix of the molecular bridge. After introducing 
the ferromagnetic Anderson-Holstein model in Sec. H, 
we discuss the necessary modifications in the quantum 
kinetic equation due to the spin-polarized leads, as op¬ 
posed to the case of non-magnetic leads. Most impor¬ 
tantly, we find an exchange magnetic field which acts on 
the dot electrons and which is generated by the coupling 
of the electrons to the vibrational mode. We discuss this 
effect for realistic experimental parameters in Sec. HI. 
We calculate the relevant spin observables and discuss 
spin precession and accumulation also in terms of spe¬ 
cific elements of the reduced density matrix in Sec. IV. 
Transport observables (such as the charge current) for a 
nonequilibrium situation as well as for different polariza¬ 
tion angles of the leads are presented in Sec. V. We show 
that the phonon-induced exchange magnetic field gen¬ 
erates a current blockade which is weaker as compared 
to the pure spin-valve effect. To complete the physical 
picture, we also provide the results for the general case 
when both the Coulomb and the electron-phonon cou¬ 
pling are present. We underpin these numerical findings 
in terms of an attractive effective Anderson model for 
the polarons. We conclude and summarize our findings 
in Sec. VH. 
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II. MODEL AND DIAGRAMMATIC 
TECHNIQUE 

In order to model the Anderson-Holstein spin-valve ge¬ 
ometry, we use the Hamiltonian 

^ + Ht + Hieads- ( 1 ) 

The central region, see Fig. 1(a), is denoted as ^dot\ or, 
equivalently, as ^molecule^ in what follows and may be 
given as, e.g., a single electronic level of a quantum dot 
or a molecule. Its Hamiltonian reads {h=l) 

Hmol cjd^da- + U -f + \x • 

cr 

( 2 ) 

The fermionic operators d(j{dy) annihilate (create) an 
electron with spin cr = on the dot. Moreover, the 
electron number operator on the dot is = d\dcr. The 
single particle energy is given by Ccr and may be con¬ 
trolled by an applied gate voltage and/or an external 
magnetic field. In the absence of an additional exter¬ 
nal magnetic field, we have that = cq. Onsite 

Coulomb interaction is taken into account by the param¬ 
eter U. Furthermore, the vibrational mode is character¬ 
ized by the bosonic operators 6 / 6 ^, and its frequency is 
denoted by Q. The (spin-independent) electron-phonon 
coupling strength is provided by A. The oscillator dis¬ 
placement operator is denoted by x, and the harmonic 
eigenstates are labeled by the quantum numbers \m). In 
the absence of the coupling to the phonon, A = 0 , the 
electronic eigenstates {| 0 ), |t), |i), \d)} correspond to the 
empty, the two singly and the doubly occupied state of 
the quantum dot. 

Spin-dependent transport features emerge due to fer¬ 
romagnetic leads. We allow for different hybridizations 
for the up- and the down-spin electrons between the leads 
and the dot and describe them by the Hamiltonian 

Hleads — E (^kr A^a) ‘ (^) 

kCKT 

Electrons in lead a = L/R with energy Ckr are cre¬ 
ated (annihilated) by the fermionic operators cj^cK^(ckar)5 
where r = ± labels the majority/minority spin species. 
The quantization axis of the left and right lead, respec¬ 
tively, is defined by the respective magnetization direc¬ 
tions rvL and xvr. Both enclose an angle (pa with the x- 
axis, see Fig. 1(b). We denote by the electro-chemical 
potential of the lead a. Ferromagnetism is included by 
assuming different densities of states pr=± at the Fermi 
energy for majority/minority spin electrons in the respec¬ 
tive lead. The polarization of the ferromagnetic lead a 
is defined by Pa = • For nonmagnetic leads, 

Pa = O 5 whereas Pa = I corresponds to a fully polarized 
lead hosting majority spins only. The dot and the leads 
are tunnel coupled with an amplitude leading to the 
Hamiltonian^^’^^ 

Ht=Y 1 + H.c] , (4) 

kacrr 


where the unitary operator describes the rotation 
of the spin coordinate system of the on-dot electron to 
the common reference frame. As in Ref. 27 and 28, we 
choose the quantization axis for electron spins on the 
dot along the z-axis of the system that is spanned by 
ex = (nfl + nL)/|nij + hlI, By = {ur - nL)/|nfl - hlI, 
and Bz = Bx X By, see also Fig. 1(b). In passing, we note 
that for spinors which are rotated by an angle 0 about 
an axis parallel to the direction of the unit vector h, the 
following transformation holds 

gh^/ 2 )(n-a-) _ cos( 6 >/ 2 ) + ih • CTsin( 6 >/ 2 ), ( 5 ) 


where cr = {cFx^cFy^cFz) is the vector of the Pauli matri¬ 
ces. Accordingly, the coordinate system is rotated by the 
opposite angles, such that (see Fig. 1(b)) the tunneling 
Hamiltonian assume the form 


Ht = ^ ta 

ka 





+ i7.c. 


(6) 


Here, we have used two rotations in series about 
the z- and ^-axis, respectively, such that Uro- = 
g 2 ((/)a/ 2 )cr^g 2 ( 7 r/ 2 )cr^^ ^ fimtc difference in the elec¬ 

trochemical potentials between the left and right lead, 
eV = Pl — Pr^ a net charge current will flow and the 
system is in general in a nonequilibrium situation. For 
Pa 7 ^ O 5 the hybridization between the dot and the lead 
a becomes in general spin dependent and we may intro¬ 
duce the tunneling coupling strength Ta± = pa±- 

Hence, electronic states with spin a = t/i acquire a finite 
line width ^oLal‘^ due to their hybridization 

with the reservoirs. 

Before we proceed with the calculation of the nonequi¬ 
librium density matrix, we apply a polaron transfor¬ 
mation, such that the electrons and the phonon de¬ 
couple according to H' = Here, S = 

Ea By this, the eigenenergies of the 

bare electronic system are shifted by the polaron energy, 
such that e^r ^ e^- —A^/Q. Successively, also the Coulomb 
repulsion strength gets shifted: U U — 2\^/Vt. In pass¬ 
ing, we note that typical terms of the tunneling Hamilto¬ 
nian transform under S according to d^ 
and, correspondingly for the creation operators, see 
Ref. 22 for the details. 

Since we are interested in the physical properties of the 
dot in the first place, the lead and the bosonic degrees 
of freedom are traced out, thereby yielding the reduced 
density matrix of the system. In order to proceed ana¬ 
lytically, we exploit the fact that the lead and phonon 
Hamiltonians are quadratic in the fermionic and bosonic 
operators. By performing the trace, we implicitly as¬ 
sume that both of these subsystems remain in thermal 
equilibrium, even in the presence of a finite bias volt¬ 
age. Otherwise, dynamic nonequilibrium occupations of 
the phonon states and/or the lead electronic states have 
to be explicitly taken into account as well before tracing 
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over them. The reduced density matrix of the system 
in general has six non-vanishing entries and includes the 
coherences due to the coherent spin evolution driven by 
the ferromagnetic leads. By following the steps described 
in Refs. 37 and 38, we may obtain the kinetic equation 
for the elements of the reduced density matrix 


^(>g(<) = “ Kx,Pi‘) (‘) + 


E/ (7) 

-f —oo 


X1X2 


with the irreducible self energy kernels (t, t'). 

They include transitions of the system between even 
and odd parity states that are induced by the tunnel¬ 
ing Hamiltonian. The internal coherences are covered 
by the first term and are generated by the molecu¬ 
lar Hamiltonian They are given by the Hamil¬ 
tonian matrix elements of Hmoh with y, y' G 

[{|0), It), It), |d)} (8) |m)], of the dot-plus-phonon-system 
in the absence of the coupling to the leads. We note 
that for a spin-independent electron-phonon interaction 
A, a contribution to zeroth order in the hybridization Tq, 
arises, which is diagonal in the basis of the phonon. 

To determine the irreducible self energy kernels, we 
employ the real-time diagrammatic technique^^’^^ which 
gives a systematic expansion of the kinetic equation in 
orders of The first non-vanishing order includes all 
terms which change the number of electrons on the dot 
by one. For practical purpose, we shall include the max¬ 
imal number n of states in the bosonic Hilbert space and 
verify convergence with respect to increasing n. Then, all 
4n diagonal elements of the density matrix, are in general 
nonzero and represent the occupations of the respective 
polaron states of the dot. In addition, the overlap be¬ 
tween the singly occupied states with spin-up and -down 
is finite in the presence of ferromagnetic leads^^’^^. In the 
presence of the vibrational degree of freedom, the appear¬ 
ing Fermi functions are dressed by the Franck-Condon 
matrix elements, see Refs. 37 and 38. They simplify if 
we assume the vibrational mode to be always at thermal 

equilibrium^^d6,22,23,36^ Then, we find 



duj 

27r 


/^(w + e)C^(w), 


(8) 


where Cq{uj) = 27r Pn(A/f2)(5(ci; — nft) (see Ref. 

23) is the Fourier transform of the correlator 

c+{t) = ^e-(Vn)(6te^-*-6e--*)e-(Vn)(6t-6)^ ^ 

Moreover, we have the detailed balance relation = 

For completeness, we comment that the 
Pn{x) have an expansion in terms of modified Bessel 
functions^^. For practical reasons, infinite sums are trun¬ 
cated once convergence is achieved for the observables of 
interest. Typically n = 20 is sufficient for large A/F 10 
as studied in the following. 


Charge current and spin expectation values are the in¬ 
teresting observables, which are also accessible in exper¬ 
imental setups. In the following, we focus on the sym¬ 
metric case, Tl = F^^, and //l = = eR/2 as well as 

Pl = Pr = P and (I)l = — = 0/2 (a generalization is 

straightforward). We present results for the differential 
conductance dl/dV below, derived from the symmetrized 
tunneling current I = {II — Ir)I2^ with 

Ia{t) — 5^(^kacr^kacr)- (10) 

kacT 


We note that the current operator has a similar struc¬ 
ture as the tunneling Hamiltonian. Within the diagram¬ 
matic formalism, we have to replace a tunneling vertex 
by a current vertex and use the corresponding transition 
rates t'), for details see Ref. 38. Information 

about the spin state of the system is obtained from the 
reduced density matrix , which is related to the vector 
of the average spin <j = (cr^,, cr^, cr^). We then have for 
the averages of the spin projections 


(o-rr) = 


t , t 
P[+ P\ 


{<^y) = 


t 4- 

.Pi-Pi 


> {<^z) = 


t 4- 

p\-pI 


( 11 ) 


III. EXCHANGE MAGNETIC FIELD 

Ferromagnetic leads induce a coherent spin dynamics 
on the dot although the quantum dot itself is nonmag¬ 
netic. In lowest order in the lead-dot tunneling coupling, 
a contribution to an effective exchange magnetic field is 
generated from each of the two leads. It points in the 
respective direction of the lead magnetization. We fol¬ 
low the derivation in Refs. 27 and 28 and explore the 
fact that the charge as well as the spin degrees of free¬ 
dom obey rate equations. The dynamics of the spin can 
be formulated in terms of a time evolution equation for 
the spin vector cr{t) which describes the dynamics of the 
spin of the confined electron. This vector equation can 
be decomposed into three parts according to 

da / da\ f da\ f da\ 

\dt Kdt \dt 

where accumulation, relaxation and rotation of the spin 
are identified by the three terms. Correspondingly, the 
dynamical equation of motion for the charge degrees of 
freedom can be decomposed in a similar manner. As 
opposed to Refs. 27 and 28, all energies or differences 
of energies between singly and doubly occupied states 
have to be replaced by the corresponding energies of the 
respective polaron states, see Sec. H. A closer look to the 
third term in Eq. (12) yields the representation 

= a X (Bl + Br) (13) 

rot 

in terms of an effective exchange magnetic field with the 
contributions It is well known that for standard 
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FIG. 2. Effective exchange magnetic field in units of F in¬ 
duced by the electron-phonon and Coulomb interaction. The 
color scale shows the sum Bex — Bl + Br. Panel (a) shows 
the result for pure Coulomb interaction U = SOksT without 
electron-phonon coupling A = 0. In (b), the electron-phonon 
coupling is set to A = lOksT, Q = bksT, while U = 0. The fi¬ 
nite A broadens the resonances, shifts it to larger values of the 
gate voltage and reverses its direction as compared to panel 
(a). In (c), both interactions are present. The remaining 
parameters are ksT = F, (j) — 7^12 and p = 0.9. 


spin-valves without a vibrational mode, the Coulomb in¬ 
teraction U between electrons can generate a finite ex¬ 
change magnetic field (see Eq. (14) below for A = 0). In 
presence of the electron-phonon coupling, the vibrational 
coupling generates an effective attractive Coulomb inter¬ 
action. In turn, the resulting exchange magnetic field 
depends on A as well as on U. Overall, we find the ana¬ 
lytic expression for the exchange magnetic field 


the direction opposite to that of the lead magnetization. 
Moreover, it vanishes exactly when U = 2A^/fl. The ex¬ 
change magnetic field is shown in Fig. 2 as a function 
of the bias and the gate voltages in a color scale plot in 
units of F for three different cases. In Fig. 2(a), only 
pure Coulomb interaction is present while the electron- 
phonon coupling is switched off {U = SOksT^X = 0). 
In turn, in Fig. 2 (b), only pure electron-phonon in¬ 
teraction is present while the Coulomb term is absent 
{U = 0, X = lOksT). Here, we have chosen the frequency 
of the oscillator as = bksT. Finally, in Fig. 2 (c), both 
interactions are present {U = SOksT^X = lOksT). In 
the presence of the Coulomb interaction alone and for a 
fixed bias voltage, the exchange field shows sharp reso¬ 
nances. In between these resonances, the field decreases, 
goes through zero, reverses its sign, and increases again 
to the adjacent resonance in Fig. 2(a). For pure electron- 
phonon coupling and in absence of a Coulomb interaction 
U = 0 (see Fig. 2(b)), the analogy to an effective attrac¬ 
tive Coulomb interaction appears in two respects: First, 
there is a finite exchange field induced by the phonon on 
the dot. Second, as opposed to (a) it has a reversed sign. 
As in (a), we find two resonances, which are shifted to 
higher values of eo here. This is a general feature which 
is due to the opposite sign of the effective Coulomb inter¬ 
action constant. In panel (c), the parameters are chosen 
such that A^/f2 « U. The exchange magnetic field is still 
visible, although being smaller in amplitude as opposed 
to (a). The width of the resonance is broader than in 
(a) which reflects the dominance of A. For those param¬ 
eters, we have that X‘^/Q > U{2 and the tendency is as 
in Fig. 2(b). In contrast, the choice 17/2 > A^/f2 would 
result in a reversed direction of the effective exchange 
field again. In what follows, we discuss how the results 
are modified in comparison to the case U 7 ^ 0 , A = 0 by 
a finite electron-phonon coupling. 


IV. SPIN PRECESSION AND SPIN 
ACCUMULATION 


B 


a 


Pa^a 

TT 




■An 





Ha, 


(14) 

which acts from contact a on the spin of the 
electron confined on the dot. Here, Aq,(x) = 

Reff^ [1/2 + ^/^(x —/iQ,)/(27r)], with the Digamma func¬ 
tion 4^^^. It is important to realize that even in the 
absence of the Coulomb interaction, i.e., for U = 0 , a 
finite exchange field Bq, in the direction of the mag¬ 
netization of the lead a exists when an electron-phonon 
coupling is present, A 7 ^ 0 . Consequently, the spin of 
a confined electron, which is subject to electron-phonon 
coupling, will process in the effective exchange magnetic 
field generated when A 7 ^ 0. Eq. (14) includes previous 
results in the case of A = 0. Furthermore, it is inter¬ 
esting to realize that for Coulomb interaction strengths 
U < 2A^/0, the exchange magnetic field points into 


The immediate consequence of the interaction-induced 
exchange magnetic field is a spin accumulation and spin 
precession of single electron spins on the dot. Whenever 
the dot is doubly occupied, these effects are expected to 
be less pronounced, since the spin singlet state is not 
sensitive to inhomogeneous magnetic fields. Spin mea¬ 
surements have been performed in order to determine the 
strength of the exchange magnetic field, e.g., in carbon 
nanotubes contacted by ferromagnetic leads"^^. In the 
following, we address the spin expectation values along 
the three spatial directions with the coordinate system 
being defined as indicated in Fig. 1(b). We show results 
for the angle (p = irj2 between the left and right magneti¬ 
zation directions. This corresponds to the configuration 
ul a iii^. From Refs. 27 and 28 it is known that for small 
bias voltages, a finite spin accumulation occurs along the 
direction np — ur = e^. It also occurs in the absence of 
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FIG. 3. Interaction-induced spin accumulation illustrated by 
the spin expectation values {(Tx) (left column), {(Ty) (middle 
column) and {<7z) (right column) for varying bias and gate 
voltages and for U = 30k bT. The upper row shows the re¬ 
sults for pure Coulomb interaction without electron-phonon 
coupling A = 0. The lower row refers to the case with an ad¬ 
ditional electron-phonon interaction A = lOksT^Q, = ^ksT. 
The remaining parameters are as in Fig. 2. 


the exchange field. It is important to realize that a finite 
(o'x) 7 ^ 0 and/or {az) 7 ^ 0 is a clear signature of a finite 
exchange field generated by Coulomb or electron-phonon 
interactions. 

The induced spin accumulation is illustrated in Fig. 3 
by the spin expectation values {(7x)^{(7y) and {(Tz) for 
varying bias and gate voltages. The upper row shows 
the results for the case of pure Coulomb interaction, 
U = 30kBT and A = 0. For 0 = 7 r/ 2 , spin is accumulated 
in the direction 11 ^ — 117 ^ due to the conservation of the to¬ 
tal angular momentum. This maximal spin-valve effect 
competes with the spin precession around the induced 
exchange field which points in the direction -f rii^. 
The latter induces finite x— and ^-components of the 
spin. Since two addition energies are present, also two 
distinct resonance lines appear in the upper left panel. 
The same substructure also emerges in the upper middle 
panel. In this case, the transitions on the dot between 
the states |t) and |^) and the empty or doubly occupied 
state are non-degenerate due to the finite spin accumula¬ 
tion. In the lower row of Fig. 3, we depict the results for 
U = lOksT and A = lOksT which illustrate the modifi¬ 
cation of the spin accumulation and the spin precession 
due to a finite electron-phonon coupling A. They are in 
accordance with Fig. 2(c), where a finite exchange field is 
shown. Although the spin expectation values are about 
one order of magnitude smaller than the exchange field, 
the qualitative physical behavior is the same. The elec¬ 
trons couple to the phonon mode differently depending on 
their spin due to the lead polarization. By this, an effec¬ 
tive spin-phonon coupling emerges. On the other hand, 


FIG. 4. Same as in Fig. 3, but as a function of the enclosed 
angle 0 for a single-particle energy eo = 30k bT and for U = 
30k bT. The upper row shows the results for A = 0, while the 
lower row shows the case with an additional electron-phonon 
interaction A = lO/c^T, Q = 

the electron-phonon coupling compensates the repulsive 
Coulomb interaction and allows the dot to be occupied 
by two electrons. 

A. Angular dependence 

The choice of two independent magnetizations in the 
left and right lead, characterized by the relative angle 
generates two non-collinear magnetic fields for the spin 
of an electron on the dot. When the dot is occupied by a 
single electron, the spin state is sensitive to these fields. 
Then, also the dependence of the spin accumulation on 
(j) is relevant. In Fig. 4, we show the spin expectation 
values as a function of the bias voltage and the enclosed 
angle (j). We have fixed the gate voltage to eo = 30k bT 
and the remaining parameters are as in Fig. 3. Again, 
the three upper panels refer to pure repulsive Coulomb 
interaction with U > 0. The device is largely singly oc¬ 
cupied and we find a clear precession of the spin with 
finite x- and ^-components of the spin operator. In con¬ 
trast to the spin accumulation, the spin precession van¬ 
ishes for the collinear configurations, when 0 = 0 , Ttt. 
As the spin-dependent effects are absent for zero bias 
voltage, we find that increasing the bias voltage favors 
a specific spin direction, either up or down. This is ac¬ 
companied by the precession of the spin with the same 
or the opposite components of the spin vector. This is 
depicted in the left and right panels of the upper row in 
Fig. 4. The results for finite A are shown in the lower 
row. The spin accumulation is sensitive to the change of 
the angle as well, similar as the behavior discussed above. 
However, the characteristics of the precession is different, 
since it is determined by the value of the exchange mag- 
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FIG. 5. Differential conductance of the Anderson-Holstein 
spin-valve as a function of the bias and the gate voltage in 
units of Gtunn = (e^//i)r/ ( 27 rA;sF). In (a), the noninteracting 
case is shown: U = A = 0, while (b) refers to U = 30k bT, X = 
0. In panel (c), we set U = 0, A = lOksT, and in (d), we have 
chosen U = 30/csF, A = lOksT^Q = bksT. All calculations 
have been performed with Q, = ^ksT, p = 0.9 and 0 = 7r/2. 

netic field. We note that always a positive x-component 
of the spin emerges, whereas the ^-component may re¬ 
verse its sign, when either the angle or the bias voltage 
are changed from positive to negative values. Since the 
Franck-Condon blockade is well developed in the central 
areas of the panels, spins of any kind are blocked to enter 
the dot in this regime. Hence, neither spin accumulation 
nor spin precession occurs here. 


V. TRANSPORT SPECTROSCOPY 

Measurements of the current give precise insight into 
the spectrum of the interacting quantum dot. In this sec¬ 
tion, we present results for the differential conductance 
dl/dV^ with the current defined in Eq. (10). The differ¬ 
ential conductance as function of eV and cq is shown in 
Fig. 5 (a) for [/ = 0,A = 0,F = ksT, and 0 = 7r/2. A 
single conductance line is visible since the system pos¬ 
sesses only one addition energy Eadd = ^\d) ~ “ 

^lt/1) ~ ^|o) = ^ 0 - The finite conductance is accompanied 
by the inset of the spin accumulation on the dot, see 
the discussion related to Fig. 3. In Fig. 5 (b), we depict 
the result for U = dO/c^T, A = 0. Without the ferromag¬ 
netic leads, two addition energies exist in this case, cq and 
eo + t/. Then, for cq < lOF, the dot is in the doubly occu¬ 
pied state, whereas for increasing gate voltages, the singly 
occupied state is favored. This gives rise to the diamond 
in the center region of Fig. 5 (b). For large gate voltages, 
^0 ^ 30F, the dot is empty again. The line widths are 
determined by the temperature, and we find pronounced 
satellite lines with negative values. This negative dif- 


FIG. 6. Differential conductance in the crossover regime (see 
text) for U = SOksT.X = Q = lOksT. The remaining pa¬ 
rameters are as in Fig. 5. 


ferential conductance originates from a partial blocking 
of spins, a competition between spin accumulation and 
interaction-induced spin precession. Spin accumulation 
hinders electrons, which do not have the right spin pro¬ 
jection, to enter the drain, whereas spin precession coun¬ 
teracts to lift the blockade. The specific positions of the 
lines are obtained from the shape of the induced exchange 
field, see the discussion in Sec. III. With the lead magne¬ 
tizations chosen such that 0 = 7r/2, a maximal spin-valve 
effect is present. The charge degeneracy points are lo¬ 
cated at Teo = —17/2, where a maximal current flows. In 
panel (c) of Fig. 5, we show a typical Franck-Condon dif¬ 
ferential conductance. There is a large range of bias and 
gate voltages {\eV\ < AOksT^SOksT < |eo| < bO/c^T), 
where transport is largely blocked. Once the bias volt¬ 
age provides enough energy, which is of the order of the 
polaron energy A^/f2, electrons are able to tunnel through 
the dot and a current can flow. We see equidistant lines 
in the spectrum that are associated with the energy dif¬ 
ferences between neighboring states. The regular ap¬ 
pearance could be explained in terms of Franck-Condon 
parabolas that allow or forbid transitions from one vibra¬ 
tional state to another, when at the same time the num¬ 
ber of electrons is changed^^’^^’^^. Due to the polaronic 
energy shift cq ^ cq + charge state degeneracies 

are shifted differently as compared to the situation shown 
in Fig. 5(b). In particular, the shift occurs to positive 
values of cq. The effective attractive Coulomb interac¬ 
tion tends to put the dot in the doubly occupied state. 
Clearly, features of a negative differential conductance as 
in Fig. 5 (b) are absent then. In Fig. 5 (d), we depict re¬ 
sults for the scenario that both interactions are finite, i.e., 
U = SOksT and A = lOksT. For this particular choice of 
parameters, the polaron energy approximately equals the 
Coulomb repulsion and we obtain a competition between 
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FIG. 7. Angular dependence of the differential conductance 
for U = SOksT without electron-phonon coupling A = 0 in 
(a), and with A = lOksT^Q, = ^ksT in (b). The remaining 
parameters are as in Fig. 5. 
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the effects of the regimes described in panel (b) and (c). 
The parameters are A = O = 7/ = SOksT. De¬ 

creasing the polaron energy further results in a crossover 
to the behavior shown in Fig. 5(b). There, the nega¬ 
tive differential conductance is fully developed, whereas 
the Franck-Condon steps vanish successively. In Fig. 6, 
we depict the results for the crossover regime, where we 
find a superposition of modulated negative differential 
conductance lines (blue) and Franck-Condon steps in the 
dl/dV-cmve. 


FIG. 8. (a) Gomparison of the differential conductance for 
equilibrated and non-equilibrated phonon modes. Within a 
voltage range \eV\ < 20F we do not find significant devia¬ 
tions from the calculation where the phonon is assumed to 
be in equilibrium, (b) Same as in (a) but for the nontrivial 
spin projection along the z-axis. We note that (cr^) ^ 0 re¬ 
flects the finite exchange magnetic field from the leads. Other 
parameters are choses as A = 2.5r, ksT = F, 0 = — 

5F,e = -40F. 


A. Angular dependence 

The transport spectrum dl/dV as a function of the 
noncollinearity of the lead magnetizations (j) is shown for 
a fixed gate voltage cq = SO/cbT in Fig. 7. For a van¬ 
ishing electron-phonon coupling, dl/dV{(j)^ eV) for A = 0 
is shown in Fig. 7(a). The dot is essentially unoccupied. 
The symmetry of the figure is due to the symmetry of 
the underlying setup, and we recall that we calculate the 
symmetrized current. For small bias voltages, a positive 
conductance peak appears for varying (j). Along 0, we 
recover the known suppression of the conductance away 
from the maximum according to sin^ 0, see Refs. 27 and 
28. With increasing bias voltage, the dot is more pop¬ 
ulated. Depending on the noncollinearity, spin is accu¬ 
mulated on the dot and processes. Transport is then 
reduced, and hence, a clear negative differential conduc¬ 
tance feature occurs. For even larger bias voltages, trans¬ 
port is independent of (j). In Fig. 7 (b), we show the 
results for U = SOksT and A = 10k bT. Again, conduc¬ 
tance lines are symmetric with respect to 0 = 0, eR = 0, 
and the strong Franck-Condon blockade is independent 
of the enclosed angle. Since the electron-phonon coupling 


favors a doubly occupied dot state, a spin singlet state 
is formed. This singlet is independent of any magnetic 
field and thus is not influenced by changes of the angle of 
the lead magnetizations. As a function of 0, the respec¬ 
tive lines show a maximum for 0 = 0, and the reduction 
follows again a ^ sin^ 0 functional form. 


VI. NONEQUILIBRIUM PHONON MODE 

In the discussion above, we have assumed that the 
phonon mode is always kept at thermal equilibrium. In 
this section, we address the situation when this no longer 
holds and the phonon mode is treated fully nonadiabat- 
ically. Starting from the Hamiltonian in Eq. (2), we 
have to truncate the phonon Hilbert space at a cer¬ 
tain number K of energy eigenstates. When calculat¬ 
ing the kernels within the real-time diagrammatic ap¬ 
proach, we explicitly keep the dependence on these basis 
states. Clearly, the truncation of phonon absorption and 
emission processes becomes increasingly cumbersome for 
large electron-phonon coupling. Especially in the con¬ 
text of cooling or heating the environment, see Ref. 34, 
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the full density matrix of the electron and phonon de¬ 
grees of freedom is of importance in order to monitor the 
effective temperature which refers to the mean energy 
contained in the vibrational mode. In order to include 
nonequilibrium effects in the phonon subspace, we trace 
over the bosonic degrees of freedom only when we com¬ 
pute electronic observables, such as the current or the 
spin expectation values. Otherwise, we keep the {AK)^ 
elements of the density matrix of the dot. We show in 
Fig. 8 (a) the differential conductance as a function of the 
bias voltage for A = 2.5r,f^ = bT^kBT = F, 1/ = SOF, 
and e = —dOF. The solid line refers to the results with the 
assumption of a thermally equilibrated phonon, whereas 
the dashed lines mark the results for different numbers 
of phonon states that are taken into account. The red 
(blue) line is for K = 9{K = 10) phonon states. We 
find a good agreement for small to intermediate bias 
voltages 0 < \eV\ < SOF. Small side peaks develop at 
larger bias voltages, which are associated to the influ¬ 
ence of the nonequilibrium phonon distribution. How¬ 
ever, the key features, such as the appearance of the neg¬ 
ative differential conductance and the Lorentzian peaks 
at \eV\ = U/2 are reproduced by the equilib¬ 

rium phonon distribution as well. The height of the 
peaks differ to some extent between the equilibrium and 
the nonequilibrium phonon model. For larger values of 
A > SF, it is difficult to obtain numerically converged 
results and the calculation in the Franck-Condon regime 
become increasingly cumbersome. In Fig. 8 (b), we com¬ 
pare, as an example, the ^-component of the spin for both 
models. Again, the qualitative form of the spin expec¬ 
tation value as a function of the bias voltage agrees for 
both types of calculations, apart from small deviations 
and different peak heights. 


VII. CONCLUSIONS 

The interplay of ferromagnetism. Coulomb and 
electron-phonon coupling in quantum dots or molecular 
transistors gives rise to interesting spintronic effects un¬ 
der nonequilibrium conditions. In the presence of ferro¬ 
magnetic leads, the nonmagnetic molecular bridge sup¬ 
ports a coherent time evolution of the spin between the 
It) and It) states. When the tunnel coupling is weak, 
we obtain the stationary density matrix in the presence 


of a finite bias voltage by means of the real-time di¬ 
agrammatic technique. Thereby, we include all quan¬ 
tum coherences that govern the coherent evolution of the 
confined electron spin. As a central result, we find a 
phonon-induced exchange magnetic field which acts on 
the electron spins confined in the intrinsically nonmag¬ 
netic dot. The electron-phonon coupling induces an ef¬ 
fective attractive Coulomb interaction, which, together 
with the spin-valve effect, induces a precession of the 
confined electron spins in the dot. Interestingly enough, 
we find that the induced exchange magnetic field points 
in the opposite direction when A 7 ^ 0 , as compared to 
1/ 7 ^ 0. The phonon-induced exchange field shows much 
broader resonances than the Coulomb-induced exchange 
field. This behavior is related to the small energy spac- 
ings between the adjacent energy levels in the polaron 
picture. However, once a single spin is confined on the 
dot, it starts to process in the effective exchange field. 
When A is finite, Franck-Condon physics emerges. In 
particular, the Franck-Condon blockade give rise to an 
empty dot, and counteracts the spin-valve physics. Nev¬ 
ertheless, if both interactions are of the same order of 
magnitude, spin-valve physics survives, and spin accu¬ 
mulation together with a spin precession is observable. 
These effects are to some extent weaker as compared to 
the standard pure spin-valve scenario. Also the sym¬ 
metrized tunneling current and the differential conduc¬ 
tance are accessible. In addition to the known results for 
the conventional quantum dot spin-valve and the estab¬ 
lished Franck-Condon sidebands, we propose a combined 
d//dU-diagram, where the negative differential proper¬ 
ties of the quantum dot spin valve are modulated by 
the Franck-Condon sidebands. These theoretically pro¬ 
posed effects might be observable in ultraclean carbon 
nanotubes or gated single molecule setups, for instance, 
by magnetic scanning tunneling microscope experiments 
in the near future. 
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